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1. SCOPE 


This specification establishes the final design of modifications 
to the display processor within the EOD-LARSYS system, as speci- 
fied in IDSD category 1 task agreement, titled "LARSYS Modifica- 
tion/Fisher F“distribution Thresholding". 




2. APPLICABLE DOCUMENTS 


The following documents form a part of this specification: 


o 

0 

0 


User Documentation, EOD-LARSYS LEC-3984 Rev. 

Task Agreement Titled: LARSYS Modification/Fishcr F- 

distribution Thresholding 


TIRF 76-0074 
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3. SYSTEM DESCRIPTION 


3.1 HARDWARE DESCRIPTION 

' N/A 

3.2 SOFTWARE DESCRIPTION 


The purpose of this modification was to develop a fifth 
thresholding option for the DISPLAY processor of the EOD-LARSYS 
system. The new option uses program generated Fisher F'-distri- 
bution values as the threshold to determine assignment of pixels 
to the null (thresholded) class. 

The other thresholding options provided by the DISPLAY 
processor are: 

o NO THRESHOLD 

0 Thresholding with user- supplied (INPUT) threshold values 

0 Thresholding with "Empirical" distribution of quadratic 
function values (program generated) . 

0 Thresholding with CHI-Squared distribution values (program 
generated) 

The modification for the Fisher F-distribut ion option was 
incorporated into the existing program logic. An additional 
routine, FDIST, was generated to compute the Fisher F-distribution 
values. The LARSYS segmentation map (LARSAA) was modified to in- 
clude the FDIST routine in the display overlay. Record two of 
the MAPTAP file, created by CLSFYl, was modified to include the 
number of pixels in each training subclass (KEPPTS (N0SUB2) ) . The 
DISPLAY common block, DISPL, in the proc element CMRKIO was increased 
to include KEPPTS(60) . 
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Changes were made to the following routines in the DISPLAY 
processor; CLSFYl, DSPLAY, RHDIFS, SETUPS. These various 
elements will be discussed in the order of use. Appendix A 
contains a DISPLAY processor flowchart. The program listings 
are in Appendix B, in alphabetical order. 

3.2.1 DSPLAY 

3 •. 2 , 1 . 1 Linkage 

DSPLAY is called from the LARSYS system executive routine, MONTOR. 
DSPLAY calls SETUPS, DSPLYl , EMTHRS, FDIST, DSPLY2, and PRTFCT. 

I 

3. 2.1.2 Interface 

The DISPL common block is used as an interface among all the 
DISPLAY processor's routines. 

3 . 2 . 1 . 3 Inputs 

None 

3. 2. 1.4 Outputs 

Prints error message on printer - see Appendix D. 

3 . 2 . 1 . 5 Storage Requirements 
Storage used; Code = 307g, Data = 40g 

3. 2.1.6 Description 

DSPLAY is the driver for the DISPLAY Processor. It was modified 
to enable it to test for the FISFIER F-distr ibution thresholding 
option and to call subroutine FDIST to compute the threshold 


values . 


All other processint', .in DvSPLAY remains the same. 


3. 2.1.7 FLOWCHART 

See Appendix A. 

3. 2.1. 8 LISTINS 

See Appendix B. 

3.2.2 CLSFYl 

3 . 2. 2. 1 Linkage 

CLSFYl is called from CLSFY. CLSFYl calls RELERR, WRTFLD, 

THRESH, WRTMTX, MCHLSK. 

3 . 2 . 2 . 2 Interface 

N/A 

3. 2. 2.3 Inputs 

None 

3 . 2 . 2. 4 Outputs 

CLSFYl outputs a MAPTAP file. (See Appendix C for MAPTAP format) 
and outputs tabular listing on printer. 

3 . 2 . 2 . 5 Storage Requirements 

Storage used: Code = 1611g Data = 

3. 2.2.6 Description 

CLSFYl was modified to enable it to output the total number of 
training field pixels (KEPPTS) per subclass on the second record 
of the MAPTAP file. 


3, 2. 2. 7 Flowchart 
N/A 


3. 2.2.8 LISTING 

See Appendix B, 

3.2.3 SETUP3 

3 . 2. 3. 1 Linkage 

SETUP 3 is called by DSPLAY. SETUPS calls FIND, NUMBER, FSBSFL, 
REDIF3, WRTFLD. 

3. 2.3.2 Interface 

N/A 

3. 2. 3. 3 Inputs 

None 

3. 2. 3.4 Outputs 

Prints list of options selected and diagnostic messages. 

3 . 2 . 3 . 5 Storage Requirements 

Storage used: Code = 1464g Data = 

3 . 2 . 3 . 6 Description 

SETUPS initializes the reading of the classification results 
(MAPTAP) file and the input processor control card(s3. It also 
prints a list of the options that were selected and determines 
the coordinates of the rectangular areas which encompass the 
training or test fields required. SETUP3 was modified to enable it 
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to read the modified second record of the MAPTAP file, to recognize 
the FISHER option, and to print the message "Apply FISHER F-dis- 
tribution thresholds", 

3.2. 3. 7 Flowchart 

N/A 

3 . 2 . 3 . 8 Listing 

See Appendix B. 

3.2.4 REDIF3 

3, 2. 4.1 Linkage 


REDIF3 is called by SETUP3. REDIF3 calls CHIN, NXTCHR, FIND, 
FLTNUM, LAREAD. 

3.2.4. 2 Interface 

N/A 

3. 2.4 . 3 Inputs 

REDIF3 reads the supervisor (control) cards in the input cards 
and STAT data from the MAPTAP file. 

3 . 2 . 4 . 4 Outputs 

Prints diagnostic messages. 

3.2.4. 5 Storage Requirements 

Storage used: Code = 1452g Data - 
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3. 2.4.6 Description 


REDIF3 reads and analyzes the input processor control cards to 
set processor control flags. The routine was modified to recog- 
nize the FISHER option card and in this process, the logic to 
test the threshold options was changed because of the increased 
number of options. A new feature was added to limit the threshold 
option to be the last one read, in case more than one threshold 
option was requested. The FTSHER threshold option includes a 
test to assure that the number of pixels foi each subclass is 
greater than the number of channels, which is required for the 
threshold value to be computed correctly for the given subclass. 

3. 2. 4. 7 Flowchart 

N/A 

3 . 2 . 4 . 8 Listing 

See Appendix B. 

3.2.5 FDIST' 

3 . 2 . 5 . 1 Linkage 

FDIST is called by DSPLAY. FDIST calls FISHIN. 

3. 2. 5. 2 Interface 

N/A 

3 . 2 . 5 . 3 Inputs 

None 

3 . 2 . 5 . 4 Outputs 


• 33 = 5 “ 
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None 


3 . 2 . 5 > 5 


Storage Requirements 


Storage used; Code * 103g Data ■ 136g 
3, 2. 5,6 Description 

The FDIST routine computes the threshold values for the FISHER 
F-distribution option. FDIST uses the FISHIN routine from the 
UNIVAC system's STAT-PACK library to obtain the FISHER F^-distri' 
bution values. The following equations are used to arrive at 
the threshold value; 

VAR « FISHIN(F(I) ,N1,N^$1) ---value returned from FISHIN 
Parameters ; 


F(I) » confidence level , for subclass I 
N1 -No. of channels 

N2 = No. of samples - No. of channels (i.e., KEPPTS (I) -Nl) 
$1 = overflow condition address return. 


FK *= P>t:N+l)*CN-l) 


where P = No. of channels 
N = No. of samples 


TMRES = FK*VAR 


The threshold values are computed for each subclass. If the 
FISHIN routine returns an overflow condition for a particular 
subclass, the threshold value will be set to 999.999 for that 
subclass . 


3 . 2. 5 . 7 Flowchart 






N/A 


.5.8 Listing 

See Appendix B. 




4.0 OPKRATION 


The EOD-LARSYS system of processors operates on the UNIVAC 
EXEC2 computer cystem. The user is required to set-up an input 
deck of cards containing the processor control cards as described 
in the LARSYS users document and to submit the deck to the 
computer center. 

The following processor conti’ol card is required to Initiate 
the F-distribution threshold options. Changes were made to the 
DISPLAY processor operating instructions in the LARSYS user docu- 
ment to reflect the addition of this processor control card: 

OPTIONS . FISHER Compute thresholds from the FISHER 

* F-distribution, using the confidence 

levels input on the threshold control 
card. 


* If the threshold control card is input, one of the four options 
(CHI-SQUARK., FISHER, EMPIRICAL, or THRESHOLD VALUES) should 
be input also. If the option card is omitted and the THRESHOLD 
card is input, CHI-SQUARE is assumed. If more than one thresh- 
old option is input, only the last one road will be performed. 


If 


5.0 TF.vST PROCEDURE 


Make runs on each o£ the threshold options including the 
"no threshold" option to verify the operation and performance 
of the new FISHER F~distribiition threshold option, and test the 
integrity of the system. In addition the STATS, PLOT, PLT , 
OUTLINE and FILTER options are turned on for each run to make 
all runs on the same set of data. 

The outputs of the five thresholding options are compared and 
check for accuracy. The FISHER input and output is analyzed. 
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6,0 TEST VERIFICATION 


For LARSYS Modification/Fisher F-Distrlbution Thresholding 


This verification Is being conducted to insure that the delivered 
program products satisfy the requirements as originally stated by 
the requesting organization. 

c ■ / ' 7/^- 

NASA Monitor 




Requestor 



Developer 





Cognizant System Manager 




Verification Date: // /-^ / 7 
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00415 

115* 

TSTFLOII.OOFlOI IsINFI5 I 


00415 

115* 
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APPENDIX C 
MAPTAP FORMAT 


MAPTAP FORMAT 


The tape, MAPTAP, is output by the processor CLASSIFY. It contains 
the statistics actually used in classification, the training field, 
category, class, and subclass information, and the classified data. 

Each file consisting of the following types or records: 

4 - run header records 
Repeated for each I - field header record 
classified field N - data records 

^1 - end-of-field record 
1 - end-of-run record 
end-of-fiie mark 

All records are written with a nonformatted FORTRAN write state- 
ment . 


RUN Header Record 1 


WRITE (MAPTAP) (DATE (I) , 1=1 , 2) ,BMFLG,BMCOMB,BMFEAT,NOCLS2 , 
NOFLD 2 , NOSUB 2 , NOFET 2 , TOTVT 2 , NOCAT , VARSZ 2 , 
(FETVC2(I) ,x-l,NOFET2) 


FORTRAN NAME 

AND DIMENSION DESCRIPTION 

DATE (2) Date classification was performed 

BMFLG Flag indicating B-MATRIX was used 

in classification 

BMCOMB No. of linear combinations in 

B-MATRIX 




FORTRAN NAME 

AND DIMENSION DESCRIPTION 


BMFEAT 

NOCLS2 
NOFLD2 
NOSUB 2 
NOFET2 
TOTVT2 
NO CAT 
VARS Z 2 

FETVC2CNOFET2) 


No. of channels used in computing 

the B -MATRIX 

No. of classes 

No. of training fields 

No. of subclasses 

No. of channels used in classificatio 
No. of vertices in training fields 
No. of categories 

Size of covariance for each subclass 
Actual channels used in classificatio 

RUN Header Record 2 


WRITE (MAPTAP) (CATNAM(I) , 1=1 ,N0CAT1) , (CLSMTX(I) ,I=1,N0CLS2) , 
(SUBNO(I) ,I=1,N0CLS2) , (SUBDES (I) , 1=1 ,N0SUB2) , 
((FLDMTX(I,J) 1=1,4) ,J=l,NOFLD2) ,( (VERTEX (I ,J) , 

1=1,2) ,J=l,TOTVT2) , (SUBCAT (I) ,1=1, N0SUB2) , 

(CLSVC2(I) ,I=1,N0SUB2) , (KATNO (I) , 1=1 ,N0CLS2) , 
(KEPPTS(I) ,I=1,N0SUB2) 

FORTRAN NAME 

AND DIMENSION DESCRIPTION 

CATNAM (NOCAT 1) Category names (if available.) 

NOCATl = no. of categories if 
CATEGORY classifier was applied 
NOCATl = no. of classes if STANDARD 
classifier was applied. 

Class names 


0 


CLSMTX(N0CLS2) 


FORTPAN NAME 
AND DIMENSION 

SUBNOCNOCLS2) 

SUBDES(NOSUB2) 

FLDMTX(4,NOFLD2) 


VERTEX(2,TOTVT2) 

vSUBCAT(NOSUB2) 

CLSVC2 (NOSUB2) 

KATNOCNOCLS2) 

KEPPTSCNOSUB2) 

WRITE (MAPTAP ( (COVMTX (I 
(CAVEMTXCI 

FORTIUN NAME 
AND DIMENSION 

COVMTX (VARSZ 2 , NOSUB2 ) 


DESCRIPTION 


No. of subclasses in each class 

Subclass names 

Training field information 

1 - field name 

2 - Class number field belongs to 

3 - Subclass number field belongs to 

4 - No. of vertices in this field 
Vertices for all the fields 
Vertices are ordered (sample , line) , 
(sample , line) 2 » • • • (sample , line)rpg,py,j .2 

Contains the category number to 
which each corresponding subclass 
belongs 

Contains the class number to which 
each corresponding subclass belongs 
Contains the category number to 
which each class belongs 
Contains the total number of training 
field pixels in each subclass 

RUN Header Record 3 


J) ,1=1 .VARSZ2) J=1,N0SUB2) 
J) ,I=1,N0FET2) ,J=1,N0SUB2) 


DESCRIPTION 


Original or B-tranformed covariance 
matrix for each subclass 
Mean vector for each subclass 


AVEMTX(NOFET2,NOSUB2) 


APPENDIX D 
ERROR MESSAGES 




DIAGNOSTICS (Display Processor) 


1. *****DSPLAY/SETUP3 ... ERROR CONDITION ON ATTEMP TO 

POSITION MAPTAP OVER FILES. 

*****FSBSFL STATUS CODE = . .. ABORTING RUN **’*' 

The system routine FSRSFL for positioning files has 
encountered difficulties in positioning the MAPTAP to 
correct file. Error occured in SETUP3 routine for 
DISPLAY. 

User should make sure he has indicated the correct 
file number for the MAPTAP, and that the MAPTAP does 
in fact have the correct number of files. 

2. ***DSPLAY/SETUP3 CORE OVERFLOW --- EXECUTION TER- 
MINATED*** 

Subroutine SETUP3 computes the storage needed for the 
specific problem, if more is needed than available this 
diagnostic is printed. 

3. AA0650 INVALID SUPERVISOR CONTROL CARD. 

The invalid card is printed along with this message. 
Check the spelling of the KEYWORD. 

4. WRITE ON UNIT N TERMINATED ABNORMALLY DAS TAPE NOT 
CREATED. I ST AT = 


Attempt to write on DAS output tape failed. This usually 
indicates a bad tape. ISTAT is the status code returned 
from the system binary I/O routine NTRAN. Execution 
continues . 

S? 




5. END OF TAPE ENCOUNTERED ON DAS UNIT LAST LINE ^ N. 

The end-o£-tape marker was encountered on the DAS out- 
put tape. The last line written was N. Execution of 
DISPLAY continues without further attempts to write on 
DAS tape. 

6 . **«***FISHER threshold REQUESTED-NOT PERFORMED 
...NO. SAMPLES FOR SUBCLASS " NAME *’ (=N) IS LESS 
THAN OR EQUAL TO NUMBER OF CHANNELS (=M) 

The program compares the number of samples (pixels) 
for each training subclass to the number of channels. 
If the number of samples is less than or equal the 
number of channels, the threshold request is bypassed 
for the subclass. 

7. FDIST-OVERFLOW CONDITION IN FISHIN ROUTINE FOR 
SUBCLASS=XXXX. THRESHOLD SET TO 999.999 

The FISHIN routine has returned an overflow condition. 
The threshold value is set to 999.999 by the program, 
for the subclass. 




